Discrete saturation thickness and anomalous potential height of native ultrathin 

aluminum oxide tunnel barriers 



K. Gloos, P. J. Koppinen, and J. P. Pekola 
Department of Physics, University of Jyvaskyld, PO Box 35 (YFL), FIN-40351 Jyvdskyld, Finland 

(February 1, 2008) 

We have investigated planar metal - insulator - metal tunnel junctions with aluminum oxide as 
dielectricum. These oxide barriers were grown on an aluminum electrode in pure oxygen at room 
temperature till saturation. We observed discrete barrier widths separated by As ~ 0.38 nm, corre- 
sponding to the addition of one layer of oxygen atoms. The minimum thickness of so ~ 0.54 nm is 
due to a double layer of oxygen. We found a strong and systematic dependence of the barrier height 
<l>o on the thickness s like "I>o ~ 2.5 eV/s'^(nm), which nearly coincides with the kinetic electron 
energy E = /2ms'^ for which the deBroglie wavelength matches the width of the barrier. 
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Aluminum oxide (AlO^;) is one of the standard con- 
struction materials for tunneling barriers because it 
is easily and reliably being formed. It is important 
for applications like single-electron transistors and 
Coulomb-blockade thermometers Q. Aluminum oxide 
is also a promising candidate to replace Si02 based gate 
dielectrics in miniaturized electronic circuits, see for ex- 
ample Ref. 1^. Although many experiments take advan- 
tage of AlOa; tunnel barriers, few reports focus on their 
intrinsic properties. The conventional point of view on 
these properties summarizes as follows: i) Almost any 
conductance per area can be achieved by simply choos- 
ing the right oxidation pressure and time, ii) The typical 
height of the tunnel barrier is about 2 eV. But experimen- 
tal data vary from below 0.1 eV up to 8.6 eV in) 
The dielectric constant e of the aluminum oxide barrier 
is smaller than that of bulk AI2O3, which is around 4.5 - 
8.9 at around 295 K 0. Like the height, e does not seem 
to be a well-defined property of AlOa,. 

In view of the importance of these tunnel barriers we 
have investigated them in more detail. We show that the 
height of ultrathin aluminum oxide barriers is not at all 
related to the band gap of bulk AI2O3. Instead, it seems 
to be a purely geometrical effect set by the lowest ki- 
netic electron energy for which the deBroglie wavelength 
is smaller than the width s of the barrier. This interpre- 
tation is supported by resolving, for the first time, the 
layer-by-layer thickness of the saturated oxide barrier. 

The most simple picture of a tunnel junction assumes 
an s wide rectangular barrier of height $0 across which 
electrons tunnel. In Wentzel-Kramers-Brillouin (WKB) 
approximation the conductance per area at zero bias and 
at low temperatures becomes |^ 

= expl— V2to$oI- (1) 

The first order correction increases the conductance 
quadratically with bias voltage V like |^ 

9iV) = 50 (1 + VyVo') (2) 



with Vq = (4?i^/e^m)<i>o/s2. Equations |l] and || can be 
used to estimate s and $0 of the junction. In general, 
image forces act on the electrons as they move through 
the barrier from one metal electrode to the other. They 
distort the potential distribution of real junctions, as de- 
scribed by the standard Simmons model The barrier 
height and its width determine then, together with the di- 
electric constant e, the zero-bias conductance go and the 
slope parameter Vq. Zero-bias anomalies (ZBA) omitted 
so far are discussed below. 

We fabricated 99 samples on an oxidized Si substrate, 
using the standard two-angle shadow evaporation tech- 
nique to deposit the metal electrodes in an ultra-high 
vacuum (UHV) chamber at a rate of 0.3nm/s. The sam- 
ples had iV = 12 nominally identical junctions in series, 
each of A = 2 X 2 /im^ cross-section. Electrodes were 
evaporated through a 300 nm thick Si3N4 mask as shown 
in the inset of Fig. ||. A 25 fim diameter Al wire sep- 
arated the mask from the substrate. By using a metal 
wire as a spacer we avoid organic material near the sam- 
ple which could otherwise contaminate the oxide layer or 
the electrodes. The first Al layer, either pure Al or its 
alloys with Cu, Au, or Nb, was exposed to pure oxygen at 
ambient temperature (~ 295 K) to form a thin insulating 
AlOx barrier. About 5 • 10^ Pa oxygen pressure was ap- 
plied for 30 minutes in the loading chamber of the UHV 
system. Then the second electrode (Al, Cu, Nb, Au, or 
a mixture between two of those metals) was evaporated. 
The thickness of the electrodes was about 90 nm. We 
have found no systematic dependence of the final tun- 
nel conductance per area on oxygen pressure, oxidation 
time, evaporation rate, junction area, or thickness of the 
metal electrodes when the relevant parameters were var- 
ied by a factor of two. In this respect the experiments are 
quite reproducible and represent the saturation thickness 
of the barrier at 295 K in pure oxygen. 

Both the size of the area and the number of junctions 
are a compromise: The area is large enough to be reli- 
ably measured, but sufficiently small to show a resolvable 
Coulomb blockade signal at 4.2 K when the sample is in 
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the liquid helium of a transport dewar. We have chosen 
this condition to make a reliable comparison and to min- 
imize the turn-around time. The number of junctions N 
still allows sufficiently large voltages to be applied across 
each junction using conventional electronic equipment. 
But simultaneously it reduces the risk of external high 
voltage spikes destroying the junctions. 

From the sheet resistance of the shorted junctions - ob- 
tained either by omitting the oxidation process or when 
the tunnel barriers have been short-circuited by a large 
bias voltage - we estimate a residual electrical resistivity 
of the evaporated metals of 4 — 10 /zJ7cm at 4.2 K and a 
resistance ratio between 295 K and 4.2 K of around 1.5, 
indicating a good quality of the electrodes, both for the 
pure metals and the alloys. 

The differential conductance dl/dV was recorded by 
applying a small low- frequency (~ 100 Hz) ac voltage dV 
superposed on the bias voltage V to the junction and 
detecting the current modulation dl with the help of a 
lock-in amplifier. Results are always presented as the 
voltage drop per junction and the conductance per area 
A of the junctions g = {l/A)dI/dV. 

Figure |l| shows typical dI{V)/dV spectra of an Al- 
oxide-Cu and an Al-oxide-Au tunnel contact at T = 
4.2 K. Both have the same characteristic zero-bias 
anomalies. The essential difference between them is the 
strength of the ^^-dependence at larger voltages. It is 
basically this difference which indicates the thinner (Au 
contact) and the thicker (Cu contact) barrier. 

Since most of our junctions have normal electrodes (at 
T = 4.2 K), we rely on several observations to ensure 
quantum tunnelling: i) The typically ~ 15% reduction of 
the conductance on cooling from 295 K to 4.2 K |^ , which 
has been proposed recently as a good tunnelling indica- 
tor jl^. ii) The 1^^— dependence of the conductance. Hi) 
Several zero-bias anomalies (ZBA) appear at lower volt- 
ages. They can be attributed to inelastic electron-phonon 
scattering in the barrier We believe that the well- 
pronounced characteristic 110 meV anomaly, which we 
have regularly observed at thin barriers with s < 1.0 nm, 
indicates the high quality of those tunnel junctions, iv) 
Several Al-oxide-Nb junctions, prepared at the same con- 
ditions as the other samples, showed the expected su- 
perconducting anomalies, v) The known shape of the 
Coulomb blockade spectrum, see below. 

Partial Coulomb blockade at T = 4.2 K can be re- 
solved as long as the background zero-bias anomaly has 
a conductance maximum at = 0. Figure ^ c) shows 
that the spectra normalized to the smoothed background 
fit well the expected shape due to Coulomb blockade. 
This also indicates that the = 12 junctions are rather 
identical. In both cases the width of the anomaly ||] 
Vi/2 = 5.44 k^T/e derived from the fit curve corresponds 
to a nominal temperature of 4.4(3) K, in good agreement 
with the experimental 4.2 K. Because of the large con- 
tact area the Coulomb blockade anomaly is quite small. 



Nevertheless, its relative size can be used to derive the 
capacitance of the junctions C ^ . From that we estimate 
an average dielectric constant of e w 4 for the parallel- 
plate geometry with C = eegA/s. 

Surprisingly, most of the tunnel spectra are quite sym- 
metric like that in Fig. 0, even when the electrodes are 
of different metals, implying that the tunnel barriers are 
likewise symmetric. Therefore we did not apply the gen- 
eralized Simmons model for junctions with an asymmet- 
ric barrier but always symmetrized the spectra to 
obtain the average barrier heights. Plotting the differen- 
tial conductance versus V"^ reveals a linear relationship 
above the 110 meV anomaly from which Vq can be ex- 
tracted (the Al-oxide-Au samples had a very small 
contribution, making the analysis more difficult). Fitting 
the spectra above the ZBA using Simmons model and the 
average e = 4 allows us to determine s and $o- The fi- 
nal result depends only weakly on the absolute value of 
e if this is varied by ±50%, and the complete analysis 
indicates a temperature dependence of e and/or $0 [0. 

Figure ^ summarizes our main findings. First, the bar- 
rier height depends very strongly and systematically on 
the thickness s of the barrier like $0 oc 1/s^. These data 
clearly demonstrate that $0 is neither related with the 
band gap of bulk AI2O3 nor to the work function of the 
different metals. Second, data points accumulate at cer- 
tain values of the thickness, following a regular pattern 
with a As « 0.38(5) nm spacing as shown in the his- 
togram in Fig. 1^. Such a spacing could be expected for a 
homogeneously grown oxide layer, built up basically by 
adding one layer of oxygen atoms to an already existing 
one (the Al atoms are much smaller than the O atoms). 
The smallest observed thickness sq ~ 0.54(5) nm corre- 
sponds then to two oxygen layers as proposed theoreti- 
cally in Ref. |ll3|l . This minimum thickness is also consis- 
tent with recent observations of stable 0.6 nm thick AI2O3 
films grown on a Si substrate |^ and 0.59 nm AI2O3 films 
on Ta |lj] . The discreteness of the width s with reason- 
able absolute values strongly supports our way to derive 
the width, and therefore also the barrier height. 

We note here that only Al-oxide-Au junctions were of 
So type, all other samples with pure Al base electrodes, 
like Al-oxide-Al and Al-oxide-Cu, were mainly of si type, 
as shown in Fig. ^. Samples with Al alloys as base elec- 
trode, like CU0.5AI0.5, were mostly of S2 type, but a con- 
siderable fraction of them was of si and S3 type. The 
barrier of samples with a pure Al base electrode seems 
to be quite homogenous, while those of Al alloy samples 
could be inhomogeneous, resulting in an average thick- 
ness between si and S2, as seen in the histogram. Ob- 
viously, thicker barriers are more easily being formed on 
Al alloys with their degraded (irregular) surface struc- 
ture than on pure Al. On the other hand, the preference 
of the Au samples to form very thin barriers could re- 
sult from the fact that Au is difhcult to oxidize unlike Al 
or Cu. After oxidation, the top-most layer is probably 
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chemisorbed oxygen. When evaporating Al or Cu as the 
second electrode, these oxygen atoms form bonds with 
the metals, but desorb easily when Au is evaporated. 

Especially the thin oxide layers with sq and si have 
a well-pronounced 110 meV zero-bias anomaly as shown 
in Fig. |l|. This anomaly characterizes the surface (or 
longitudinal) phonons of crystalline AI2O3 (sapphire), see 
for example Ref. Thus in these samples the oxygen 
ions are fully bound in a regular crystal lattice and not 
only chemisorbed. Chemisorption may be the case for 
the thicker barriers. For them the 110 meV anomaly is 
strongly suppressed and an anomaly at around 50meV 
appears. A second reason for this suppression could be 
inhomogeneities in the barrier. This is consistent with 
observing noisier spectra at those junctions, indicating 
more defects. 

Several other experiments either contradict or support 
our findings. The supporting ones use a sample structure 
similar to ours, namely metal-oxide-metal tunnel juc- 
tions: Lau and Coleman investigated thermally and 
plasma-discharge grown aluminum oxide barriers while 
Earner and Ruggiero ||] sputter-deposited AI2O3 tunnel 
barriers. In both cases the absolute values as well as the 
thickness dependence of the barrier height ^o{s) agree 
well with our data. Especially the latter report extends 
our data set to significantly larger s. The contradict- 
ing experiments use different sample setups, metal-oxide- 
semiconductors juctions: Ballistic electron emission spec- 
troscopy on AlOx found 3.90(3) eV high barriers for 8 nm 
oxide thickness compared to our (extrapolated) 

40meV. Internal photoemission gave $0 — 3.25(8) eV for 
5—15 nm wide barriers According to Ref. there 
is almost no difference of the barrier height when AlOa, is 
replaced by Si02. But according to Ref. the numbers 
are different, namely 4.25(5) eV and 3.25(8) eV for Si02 
and AlOa;, respectively. 

The different results for the different kind of setups 
could indicate that boundary conditions matter. And 
this could lead to a tentative explanation for the well- 
defined 'cut-off energy' ^o{s) of our samples. The exper- 
imental $o(s) in Fig- § marks roughly the kinetic energy 
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at which the deBroglie wavelength of the electrons be- 
comes smaller than the width s of the barriers. The 
cut-off would then be simply induced by an electronic 
resonance phenomenon inside the barrier. 

We have tried to verify independently the potential 
height by measuring the current-voltage characteristics 
to higher voltages. Figure § shows the voltage Vt at 
which the tunnel current raises dramatically, exceeding 
~ 1 mA. The junctions are then usually destroyed, re- 
sulting in either a short or an open circuit. However, 
this breakdown voltage Vf, seems to be smaller than the 



one expected for field emission. Up to about $0 = 4eV 
the breakdown voltages are closely correlated with the 
potential height, thus supporting the magnitude of $0 
derived from the spectra at low bias voltages. The larger 
$0 above 4 eV requires additional considerations because 
of the obvious trend of Vb to saturate at around 2.5 V 
(the Al-oxide-Au samples had slightly, but systemati- 
cally, larger Vb than the Al-oxide-Cu or Al-oxide-Al sam- 
ples). Either these large barrier heights are wrong or we 
face a new effect. The latter is probably the case because 
of the huge electrical breakdown fields Eb = Vb/s, ap- 
proaching 5 GV/m at small s. Compared to that, typ- 
ical literature data for AlOx junctions have much lower 
Eb « 0.1-1 GV/m at most [|8|. A breakdown field which 
decreases with increasing gap between the electrodes is 
a well-known, yet not fully understood, phenomenon in 
vacuum high- voltage insulation [|l9|. The 2.5 V upper 
bound of Vb could be readily explained by the damage 
done due to the mechanical stress excerted by the electri- 
cal field. One could also speculate whether the extremely 
high tunnel current density ~ 10* A/m^ just below Vb in- 
duces a high enough electron density inside the barrier 
that could transform the insulating AlOa, into a conduct- 
ing metallic state. Whatever the right explanation, the 
large Eb confirms the excellent quality of our tunnel bar- 
riers with a very small number of defects that could act 
as weak spots. 

Finally, the systematic $o(s) dependence leads to a 
preferred transparency of the tunnel junctions for each 
saturation thickness which covers roughly two decades in 
conductance per area. A wider range of transparencies 
is difficult to achieve by just varying the thickness of the 
barriers. 
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FIG. 1. a) Differential conductance g = {\/A)dI/dV 
vs. voltage drop V for two typical tunnel junctions at 
T = 4.2 K. Solid hne Al-oxide-Cu, dashed hne Al-oxide-Au 
(not symmetrized original data). b) Second derivative 
g' = dg/dU of the two spectra in arbitrary units. c) 
Coulomb-blockade anomaly of the two tunnel junctions. The 
change of conductance due to Coulomb blockade has been 
normalized to the background conductance gr- Solid white 
lines are theoretical predictions. The vertical bar indicates a 
2 X 10"'* variation. From the size of the anomalies we esti- 
mate junction capacitances of about 0.22 pF and 0.18 pF and 
dielectric constants of e ~ 5.6 and 3.1 for the Cu (top) and 
the Au (bottom) sample, respectively. 
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FIG. 2. Barrier height $o versus width a of aluminum ox- 
ide tunnelling barriers. The solid line represents Eq. ^. Solid 
circles are the data of Ref. [5] for sputter-deposited AI2O3 bar- 
riers. Crosses are averaged data from Ref. [6] for thermally 
or plasma-discharge oxidized Al. The inset shows our exper- 
imental setup to evaporate the metallic electrodes through a 
thin Si3Ni4 mask onto an oxidized Si substrate (not to scale). 
Arrows mark the two different evaporation angles before and 
after oxidation. 
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FIG. 3. Histogram of the number of samples falling inside 
a 0.1 nm wide interval around a barrier thickness s. A total 
of 99 samples have been measured. Trace a represents the 
Al-oxide-Au samples, trace b belongs to (M, Al)-oxide-M and 
(M, Al)-oxide-(Al, M) samples (M = Al, Cu, Au, or Nb), and 
trace c to the Al-oxide-M samples. Trace d is the sum of all 
samples. Traces b, c, and d are vertically displaced. 
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FIG. 4. a) Breakdown voltage Vt versus barrier height "l>o 
and b) breakdown field Et versus barrier thickness s. The 
solid lines are Vt = 0.5"I>o/e and Eb = 1.25 (GV/m)/s^(nm), 
respectively. 
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